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ABSTRACT 

In this paper we applied Cloude's decomposition to imaging radar polarimetry. We 
show in detail how the decomposition results can guide the interpretation of scattering 
from vegetated areas. For multifrequency polarimetric radar measurements of a clear- 
cut area, the decomposition leads us to conclude that the vegetation is probably thin 
compared to even the C-band radar wavelength of 6 cm. For a forested area, we notice 
an increased amount of even number of reflection scattering at P-band and L-band, 
probably the result of penetration through the coniferous canopy resulting in trunk- 
ground double reflection scattering. The scattering for the forested area is still 
dominated by scattering from randomly oriented cylinders, however. It is found that 
these cylinders are thicker than in the case of clear-cut areas, leading us to conclude 
that scattering from the branches probably dominates in this case. 


1. INTRODUCTION 

In the analysis of polarimetric imaging radar data, one is often faced with the situation 
where some geophysical parameters must be inferred from an area that exhibits 
significant natural variability in the scattering properties. In such cases, the resulting 
average Stokes or covariance matrix differs significantly from that of a single scatterer. 
It has long been speculated that if a unique way could be found to decompose such an 
average Stokes or covariance matrix into a sum of matrices representing single 
scatterers, one would not only be able to more accurately interpret the scattering 
processes, but the problem of inferring geophysical parameters from the measured 
radar data would also be significantly simplified. 

One of the first examples of such a target decomposition technique was provided by 
Chandrasekhar [1] in his book on radiative transfer. Considering the case of lateral 
scattering of light by small anisotropic particles, he decomposed the total average phase 
matrix into the sum of a phase matrix that represents dipole scattering and a phase 
matrix that represents pure random scattering, or noise. As such, his target 
decomposition followed the same principle that breaks the Stokes vector of a partially 
polarized wave into the sum of a Stokes vector representing a fully polarized wave, and 
a Stokes vector representing a completely unpolarized wave (see Papas [2] for example). 
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Huynen [3], In his thesis published in 1970, also Introduced a target decomposition 
theorem In which he decomposed an average Mueller matrix Into the sum of a Mueller 
matrix for single a scatterer and a “noise” or N- target Mueller matrix. However, unlike 
Chandrasekhar’s decomposition, Huynen’s N-taiget Is not polarization independent, 
and has been shown to be just one of an Infinite set of such residue matrices, meaning 
that Huynen's decomposition is not unique. 

In 1988 Cloude 14] introduced a target decomposition based on an eigenvector 
decomposition of the target covariance matrix. This decomposition was shown to be 
unique and. In the monostatic case, breaks the average covariance matrix up into the 
weighted sum of three covariance matrices representing three different single scatterers, 
analogous to the decomposition of a Stokes vector of a partially polarized wave Into the 
sum of the Stokes vectors of two full polarized waves that are orthogonally polarized (2). 

In this paper we examine the application of Cloude’s decomposition to the analysis of 
polarimetric synthetic aperture radar (SAR) data. 

2. CLOUDE’S DECOMPOSITION 

Cloude showed that a general covariance matrix [TJ can be decomposed as follows: 

ITI = X ***** (1) 

f=l 

In (1). X i ; i = 1,2, 3,4 are the eigenvalues of the covariance matrix, k £ ;i = 1,2,3, 4 are 
the eigenvectors of [T], and kj means the cuftoint (complex conjugate transposed) of 
k ( . In the monostatic case, for reciprocal media, the covariance matrix has one zero 
eigenvalue and the decomposition results in at most three covariance matrices on the 
right-hand side of (1). 

In (1) the elements of the eigenvectors are the elements of an equivalent scattering 
matrix in the basis in which the decomposition is performed. For example, if the 
covariance matrix elements are expressed in the normal HV basis, then the four 
elements of the eigenvectors are simply the four scattering matrix elements In the linear 
basis. Also, since the covariance matrix Is, by definition, a hermitian matrix, it follows 
that the eigenvalues are all real, and the eigenvectors are orthogonal, therefore, one 
ran interpret (1) to mean that this decomposition breaks the covariance matrix into the 
weighted sum of covariance matrices from “orthogonal” scattering matrices. 

Cloude did his decomposition In the basis formed by the Pauli spin matrices, but 
pointed out that the decomposition can be done using any four complex matrices that 
satisfy the constraint of completeness and normalization. Here, we shall Interpret our 
results in the natural linear basis in which the radar measurements are performed. 
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Also useful in our discussions later is Cloude’s definition of target entropy, 

4 

H T= “X P t lQ g 4 Pi (2) 

t= 1 

where 


P< = 


A^ + A2 + A3 + A4 


(3) 


As pointed out by Cloude, the target entropy is a measure of target disorder, with 
Ht = 1 for random targets and H T - 0 for simple (single) targets. 


3. DECOMPOSITION OF SCATTERING FROM AZIMUTHALLY SYMMETRIC NATURAL 

TERRAIN 


Before getting into the radar data analysis, let us illustrate the decomposition using a 
general description of the covariance matrix for aztmuthally symmetrical natural terrain 
in the monostatic case. Borgeaud et aL [5] showed using a random medium model, 
that the average covariance matrix for azimuth ally symmetrical terrain in the 
monostatic case has the general form 
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The superscript * means complex conjugate, and all quantities are ensemble averages. 
The parameters C , 77 , £ and p all depend on the size, shape and electrical properties of 
the scatterers, as well as their statistical angular distribution. It is easily shown that 
the eigenvalues of [T1 are 
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(9) 
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The corresponding three eigenvectors are 
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(12) 
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In these expressions, we used the shorthand notation 


A = [£-l] 2 + 4|p| 2 


(15) 


We note that A Is always positive. Also note that we can write the ratio of the first 
elements of the first two eigenvectors as 


hlL = 
^21 
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which is always negative. This means that the first two eigenvectors represent 
scattering matrices that can be Interpreted in terms erf odd and even numbers of 
reflections. In the rest of this paper, we shall examine decompositions of the form 
shown above for different types of terrain using measured polarimetric radar data. 
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4. EXAMPLES OF CLOUDE DECOMPOSITION} 


4.1 Randomly oriented dielectric cylinders. 


Before analyzing measured data, let us examine a special theoretical case of uniformly 
randomly oriented dielectric cylinders. In general, the scattering matrix of a single 
dielectric cylinder oriented horizontally can be written as 


[S] = 


a 0'' 
0 b j 


(17) 


where a and b are complex numbers whose magnitudes and phases are functions of 
cylinder dielectric constant, diameter and length [61 . Generally, when the cylinder Is 
thin compared to the wavelength, |a| > |b) . To calculate the average covariance matrix 
for randomly oriented cylinders, one has to apply a rotation operator to [S] and average 
the results over all angles of orientation. Assuming a uniform distribution in angles 
about the line of sight, one can easily show that the resulting average covariance matrix 
for the monostatic case is of the form given in (4) with 


c = 

[3|a| 2 +3|b| 2 + 29t(a*b)]/8 

(18) 

p = 

[|a| 2 + |b| 2 + 6$(a*b)] 
3|a| 2 +3|b| 2 + 29i(a*b)] 

(19) 

77 — 
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(20) 

if — 

3|a| 2 + 3|b| 2 + 29t(a*b) 

C = i 

(21) 

It is easily shown using (9) — (14) that 
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The corresponding three eigenvectors are 
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From (19) and (20). it is easily shown that 

1 -H = *1 


(28) 


This means that for the case of randomly oriented cylinders, the second and third 
eigenvalues are the same. 


In the thin cylinder limit. |b| approaches zero, and we find that 
Pthin = 1/3 

Vt hin = 2/3 
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In this case, equal amounts of scattering are contributed by the matrix that resembles 
scattering by a sphere (or a trihedral comer reflector) and by the cross-polarized return, 
although a significant fraction of the total energy is also contained in the third matrix 
resembling scattering by a dihedral comer reflector. The entropy in this case is 0.95 
indicating a high degree of target disorder or randomness. 

On the other hand, in the thick cylinder limit. |b) approaches |a| and we find that 

Pthick = 1 < 32 > 

Vthick = 0 (33) 

In this case, only one eigenvalue is non-zero, and the average covariance matrix is 
identical to that of a sphere (pr a trihedral comer reflector). Hie entropy is 0, indicating 
no target rando mn ess, even though we have calculated the average covariance matrix 
for randomly oriented thick cylinders! The explanation for this result lies in the fact 
that when the cylinders are thick, the single cylinder scattering matrix becomes the 
identity matrix, which is insensitive to rotations. Therefore, even after rotations, we still 
added up only identity matrices, and no apparent target randomness is introduced. 

For cylinders that are neither thick compared to the wavelength, nor thin compared to 
the wavelength, 0 < |b / a| < 1. Figure 1 shows the behavior of the normalized 

the eigenvalues as predicted by (22)-(24). We notice, as mentioned before, that the 
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second and third eigenvalues are always the same magnitude, but they get smaller 
relative to the first eigenvalue as the ratio of the two scattering matrix elements 
approaches 1. 



Figure 1. Ratio of eigenvalues as a Junction of the ratio of b to a. The ratio of 
b to a depends on the cylinder dielectric constant and cylinder size compared 
to the wavelength. 


4.2 Measured scattering from clear-cut area 

Using the NASA/JPL AIRSAR system (7], the complete Stokes matrix [8] can be 
measured for each pixel In a scene. The actual measurements In general compare well 
with the general form of the covariance matrix assumed In (8). except that the four 
terms assumed to be zero are not exactly equal to zero. However, these four terms are 
typically much smaller than the other terms, so that to a good approximation, the 
assumption in (8) Is valid. Here, we shall set those four terms to zero. Table 1 below 
shows the measured parameters at P-band, L-band and C-band for a clear-cut area In 
the Shasta Trinity National Forest In California. 


The clear-cut area is covered with short shrub-like vegetation. We note that at all three 
frequencies the scattering is dominated by an odd number of reflections. At P-band the 
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even number of reflections and the cross-polarized returns are almost the same 
strength, and about half that of the odd numbers of reflections. This Is very similar to 
the thin randomly oriented cylinder case discussed earlier. As the frequency increases, 
the even number of reflections and the cross-polarized returns become more different, 
and also become a smaller fraction of the total scattering. This Is consistent with the 
randomly oriented cylinder case where the radius of the cylinder increases. The same 
conclusion is reached when considering the entropy. The highest value (most 
randomness) Is observed at P-band, and the randomness decreases with increasing 
frequency. As pointed out before, as the cylinder radius increases, the entropy 
decreases. Therefore, from the decomposition results we conclude that the vegetation in 
the clear-cut is randomly oriented, and that most of the scattering comes from 
vegetation that has branches that are thin compared to the three radar wavelengths. 


Table 1. Measured values for a clear-cut In the Shasta Trinity National 
Forest Also shown are the results of the Claude decomposition. 


Parameter 

P-band 

L-band 

C-band 

n 

0.5261 

0.5308 

0.4083 

C 

0.5642 

0.7580 

0.7159 

p 

0.0928+10.0582 

0.2324+10.1057 

0.3558+10.0440 


1.0260 

1.1615 

1.2437 

X2 

0.5382 

0.5964 

0.4722 

x 3 

0.5261 

0.5308 

0.4083 

Entropy 

0.95 

0.94 

0.88 


4.3 Measured scattering from a forested area. 


Table 2 below shows the measured parameters at P-band, L-band and C-band for a 
forested area In the Shasta Trinity National Forest in California. We note that at all 
three frequencies the scattering from this coniferous forest is dominated by an odd 
number of reflections. At P-band the even number of reflections contribution is 
stronger relative to the cross-polarized returns when compared to the case of scattering 
by randomly oriented cylinders. This is likely a result of the P-band signals penetrating 
through the canopy and then suffering a double reflection when scattering from the 
ground and then from the tree trun ks befo re returning to the radar. We notice that the 
imbalance of the even numbers of reflections and the cross-polarized returns decrease 
as the frequency increases, and also become a smaller fraction of the total scattering. A 
possible Interpretation of this is that as the frequency increases, the penetration 


214 









through the canopy decreases and the scattering from the randomly oriented branches 
becomes more Important. The highest value of the entropy (most randomness) is 


Table 2. Measured values for a forested area in the Shasta Trinity National 
Forest Also shown are the results of the Cloude decomposition. 


Parameter 

P-band 

L-band 

C-band 


0.3301 

0.3485 

0.2416 

c 

0.6529 

0.7122 

0.4685 

p 

0.2803+10.0167 

0.3950+10.0582 

0.3669+i0.0016 

Ai 

1.1566 

1.2805 

1.1873 

A 2 

0.4963 

0.4316 

0.2812 

a 3 

0.3301 

0.3485 

0.2416 

Entropy 

0.87 

0.84 

0.75 


observed at P-band, and the randomness decreases with increasing frequency. From 
the decomposition results we conclude that the penetration through this coniferous 
canopy cannot be neglected. Most of the scattering can still be contributed to the 
randomly oriented branches, however. The relative strengths of the different 
mechanisms, as well as the lower values of the entropy, suggest that the branches are 
thicker than those of the vegetation in the clear-cut area. 


5. DISCUSSION 

In this paper we applied Cloude’s decomposition to imaging radar polarimetry. The 
decompositions illustrated here have been implemented on a pixel-by-pixel basis to 
analyze polarimetric radar images. Space does not permit the inclusion of more 
examples, however. 

We have shown how the decomposition results can guide the interpretation of scattering 
from vegetated areas. For clear-cut areas, we concluded that the vegetation is probably 
thin compared to even the C-band radar wavelength of 6 cm. For forested areas, we 
noticed an increased amount of even number of reflection scattering at P-band and L- 
band, probably the result of penetration through the coniferous canopy resulting in 
trunk-ground double reflection scattering. The scattering for the forested area is still 
dominated by scattering from randomly oriented cylinders. It is found that these 
cylinders are thicker than in the case of clear-cut areas, leading us to conclude that 
scattering from the branches probably dominates In this case. 


215 








The quantitative analysis of scattering requires detailed analysis of the decompositions 
that will result from model predictions. Future work will concentrate on analyzing the 
decompositions of model predicted covariance matrices. The results can then be used 
to analyze polarimetric radar images. 
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